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Abstract 

Q\ , Recent results on minijet production in nuclear collisions at the RHIC and LHC en- 

ergies are reviewed. Initial conditions of the QGP at r = 0.1 fm/c, especially parton 
chemistry, thermalization and net baryon number-to-entropy ratio are discussed. Also, 
Q . contribution of minijets from a hard BFKL-pomeron ladder will be estimated. 

O: 

1 Introduction 

Particle and transverse energy production in the central rapidity region of heavy ion colli- 
sions can be treated as a combination of hard/semihard parton production and soft particle 
production. With increasing energies, the semihard QCD-processes are expected to become 
increasingly important. This is due to two reasons: firstly, already in pp(p) collisions the 
rapid rise of the total and inelastic cross sections can be explained by copious production of 
semihard partons, minijets, with transverse momenta p-r > po ~ 1...2 GeV |L|. This is also 
expected to happen in AA collisions at very high energies. Secondly, the semihard particle 
P^h! production scales as A^ 3 , so that for large nuclei the importance of semihard partons is 

increased further || [| The soft, non-perturbative, particle production in ultrarelativistic 
heavy ion collisions can be modelled e.g. through strings || ^, |7| or through a decaying 
strong background colour field [|J. 

The time scale for producing partons and transverse energy into the central rapidity 
region by semihard collisions is short, typically ~ l/po ~ 0.1 fm/c, where po ~ 2 GeV 
is the smallest transverse momentum included in the computation. The soft processes are 
completed at later stages of the collision, at r s ~ 1/Aq CD ~ 1 fm/c. If the density of partons 
produced in the hard and semihard stages of the heavy ion collision becomes high enough 
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- as will be the case - a saturation in the initial parton production can occur |9|, [K| |11 
and softer particle production will be screened. The fortunate consequence of this is that a 
larger part of parton production in the central rapidities can be computed from perturbative 
QCD (pQCD) at higher energies and the relative contribution from soft collisions with 2 
GeV becomes smaller. Typically, the expectation is that at the SPS (Pb+Pb at ^Js = 17 
AGeV), the soft component dominates, and at the LHC (Pb+Pb at ^/s = 5.5 ATeV) the 
semihard component is the dominant one. At the RHIC (Au+Au at y/s = 200 AGeV) one 
will be in the intermediate region, and both components should be taken into account. 

A lot of effort has also been devoted for building event generators |7|, [12| addressing the 
dominance of semihard processes in nuclear collisions at high energies. These have generated 
also new insight and very useful discussion during the recent years. Also recently, a promising 
novel approach to minijet production has been developed [ |l3|| . 
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I have divided this talk basically into two halves. In the first one, I will recapitulate 
the basic features of semihard parton production and review our latest results [|, [10], PI . 
The main goal of these studies is to find out the initial conditions for early QGP-formation 
at t ~ 0.1 fm/c, including the transverse energy deposited into the mid-rapidity region, 
chemical composition of the parton plasma, and, to study the possibility of a very rapid 
thermalization and estimate the initial net baryon-to-entropy ratio. It is vitally important 
to study the early formation of strongly interacting partonic matter, since the later evolution 
of the QGP, the final state global observables, and the suggested signals of the plasma will 
strongly depend on the initial conditions. The second half I will devote for discussion of an 
additional mechanism for parton and transverse energy production: minijets from a BFKL- 



ladder . Especially, I will estimate the maximum amount of transverse energy one should 



expect from the BFKL-minijets in heavy ion collisions. 



2 Initial conditions for QGP at r ^ 0.1 fm/c 

Hadronic jets originating from high quarks and gluons are clearly observed experimentally 
but when the partons have pp~ 5 GeV the jets become very difficult to distinguish |15j from 
the underlying event. In heavy ion collisions, where we expect hundreds (RHIC) or thousands 
(LHC) of minijets with p T ~ 2 GeV be produced, detection of individual minijets will be 
impossible. However, the semihard partons are expected to contribute dramatically to the 
early formation of QGP. The idea of multiple production of semihard gluons and quarks in 
pp and AA collisions is based on a picture of independent binary parton-parton collisions. 
The key quantity is the integrated jet cross section, 

(r}et(Vs,Po) = dp 2 T dy l dy 2 V] / dy 2 x 1 f i/N (x 1 ,Q)x 2 f j/N (x 2 ,Q)- I (s,t,u), (1) 
1 j p mi= J dt 

where xi j2 are the fractional momenta of the incoming partons i and j, and fiw{ x -iQ) are 
the parton distributions in iV (= p, A). The factor 2 comes from the fact that, in the lowest 
order (LO) pQCD, there are two partons produced in each semihard subcollision. In the 
eikonal models for pp collisions |IJ the ratio crj et / cr illclastic can be interpreted as the average 
number of semihard events in one inelastic collision. The results I will be quoting in the 
following 0] are obtained with the MRSH (ITH] and MRSD-' W7\ parton distributions with a 



scale choice Q = p^. More detailed formulation can be found in Refs. [|3], |TI|, and numerical 
evaluation of Eq. ([[]) in Ref. ||. 

The formula above is defined in the lowest order (LO), dajdt ~ a 2 . Often a constant 
factor K ~ 2 is used to simulate the effects of NLO terms. Studies of the NLO jet cross 



section dajidp^dy) ]T^| show that (with a scale choice Q = Pt and with a jet size R ~ 1) 



this is a reasonable approximation |18[ . Strictly speaking, however, a theoretical .fT-factor 
can only be defined for quantities where a well-defined, infrared-safe measurement function 
can be applied |19[. For .Bp-production in nuclear collisions, an acceptance window in the 
whole central rapidity unit defines such a function but for this acceptance criteria and for 
p T ~ 2 GeV the exact NLO contribution has not been computed yet. 
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The first estimate of the average number of produced semihard partons with p^ > po in 
an AA collision at a fixed impact parameter b can be obtained as || 

N AA {b, v^,p ) = 2T AA {b)a jct {^,p ), (2) 

and the average transverse energy carried by these partons as [|J 

E AA (b, y^,p ) = T AA (b)a ict (^,p )(E T ), (3) 

where T AA (b) is the nuclear overlap function || which scales T AA ~ A 4 / 3 , describing thus the 
typical scaling of hard processes in nuclear collisions. The normalization is J d 2 bT AA (b) = 
A 2 and, for large nuclei with Woods-Saxon nuclear densities, 7^(0) ~ A 2 /{nR 2 A ). The 
acceptance criteria imposed for the quantities <Tj Ct (y/s, po) and for a^ t (y/s,po)(ET) will be 
\y\ < 0.5, and the corresponding cuts will be made in y\ and y 2 . In Eqs. (||) and @ above, 
T AA (b)<jj et is the average number of semihard collisions and (E^) is the average transverse 
energy carried by the partons produced in each of these collisions. We fix po = 2 GeV, i.e. 
we describe the initial conditions at r ~ 1/po = 0.1 fm/c. The predictions for the central 
rapidity unit in Pb-Pb collisions at the RHIC and LHC energies are summarized in Tables 1. 
Also, contributions from gluon, quark and antiquark production are shown separately 

In the results given, we have neglected nuclear effects in parton distributions: fi/ A = 
Afi/ p . In reality, however, in the typical x-region x ~ 2p /v^ there are quite strong shad- 
owing corrections | ]20|] , especially for the LHC nuclear collisions. Also, the scale evolution of 
nuclear gluon shadowing was shown to be potentially important in the analysis in Ref . [21] . 



However, a re-analysis with the input from HERA at small-x |22j has to be performed before 
getting a solid quantitative prediction of the shadowing effects on minijet production. 

The rapid rise of the structure function F 2 (x, Q) at small values of x observed at HERA 



22] does not affect the bulk of the 2 GeV minijets at RHIC energies very much but obviously 
has quite dramatic consequences at the LHC energies. As demonstrated in ||, there is a 
clear enhancement of minijet production due to the new parton distributions. However, the 
more rapidly the gluon distributions rise, the more there should be nuclear shadowing due 
to the GLRMQ-fusions || Again, a more quantitative prediction depends on the 

scale evolution of nuclear gluon shadowing as well. 

Let us now have a closer look at the results in Table 1. There are four important 
observations. Firstly, the gluons clearly dominate both the initial parton and transverse 
energy production: the initial parton system is about 80 % glue. 

Secondly, the effective transverse area of the produced semihard partons is N AA n jp\. 
Comparing this with the effective nuclear transverse area, irR 2 A , we notice that 

U = NAA ^ p2 ° ~ 1 for LHC and £4 < 1 for RHIC, (4) 
txR a 

i. e. the parton system at the LHC at 0.1 fm/ c is already dense enough so that a saturation of 
parton production can take place [|, ||, |j~l"f . In this way, the scale po acquires also dynamical 
significance. At RHIC, since £4 < 1, saturation occurs at smaller values of pt (at r > 0.1 
fm/c), possibly in the region where pQCD cannot be trusted. This qualitative argumentation 

2 Results for da j '{dp-rdy) at y = can be found in Ref. JTo|. 
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Table 1: (a) The average numbers of semihard partons at r = 0.1 fm/c with \y\ < 0.5 and pr > 2 GeV 
in central Pb-Pb collisions, as given by Eq. (Q). Shadowing is not included and K = 2. The upper values 
are obtained with MRSH and the lower ones with MRSD-' parton distributions, (b) The average transverse 
energy carried by these partons, as predicted by Eq. (|^). 



is supported by a more quantitative, although still phenomenological, analysis of Ref. [ p~0f| , 
where we suggested that at sufficiently large energies (LHC) and large nuclei (A ~ 200), a 
dynamical screening mass is generated, causing a saturation in the minijet cross sections |10| 



at a perturbative scale like p ~ 2 GeV ^> Aqcd- The consequence is that the softer parton 
production is screened and its relative importance becomes smaller^] 

The third interesting observation is that the gluonic subsystem in the central rapidity unit 
Ay = 1 may thermalize very fast, at least in the LHC nuclear collisions. In the perturbatively 
produced system the (transverse) energy per gluon is E^/Np hPh = e pQ CD / n pQ CD ~ 3 GeV 
and the energy density of the system at = 0.1 fm/c is eP QCD = E^/{nR 2 A ThAy). The 
temperature T eq of an ideal (massless boson) gas in a complete thermal (= both kinetic and 
chemical) equilibrium with this energy density can be computed from e^ deal = 37r 2 /90-16T e 4 q = 
e P QCD , and we get T eq = 0.988 (1.14) GeV with the MRSH (MRSD-') densities. Especially, 
we find [(§ 

gPQCD gideal ^pQCD ^ideal 

" - ■" for LHC and -|_ > -S— for RHIC, (5) 



pQCD ideal " pQCD ideal 

Hg Ug Tig Ug 

so that at the LHC the average energy of gluons is already as in an ideal gas in thermal 
equilibrium. Only isotropization is needed, and a rapid thermalization is indeed possible. An 
instant thermalization would in turn have profound consequences on e.g. thermal dileptons, 
for which a high initial temperature plays a crucial role |25| .[] 

Note that our conclusion of the possibility of an almost instant thermalization is due to the 
small-x enhancement in the HERA gluon densities. From the energy/gluon viewpoint it also 
seems that thermalization for RHIC is going to happen somewhat later. Note however, that 
above I did not consider isotropization of the system at all. In the simplified picture presented 
here, the transit time of the colliding nuclei, r tr ~ 2R^/j, and the initial parton spread, Az ~ 
l/(xp) ~ 1/po for the partons which will be produced in the mid-rapidity, are neglected. 
Then a Bjorken-like boost-invariant picture is possible, and in the central rapidities a proper 
time t is a relevant variable. For a more thorough discussion of isotropization, a more 



detailed microscopic space-time picture has to be specified, as done in Refs. [26| (see 
also the discussion in |TT|, [TD|] ) . 



3 A similar saturation effect is also expected in the approach by McLerran et al. jl3| . 
4 On the other hand, for the thermal dileptons the trouble is the small out-of-equilibrium quark- antiquark 
component of the early parton system. 
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The fourth observation JTI]| is that initially, at r ~ 0.1 fm/c, the net baryon number 



density in the central rapidity unit is very small as compared to the gluon density but larger 
than the nuclear matter density (0.17 fm -3 ), even though the colliding nuclei are practically 
already far apart, especially at the LHC where r tr <C r^. More precisely, we estimate 

= 1 _ _ l(N q -Ng) _ j 0.25 (0.30) far 3 , for LHC 
n B -B - 3 K ng) - nR 2 AAy/po - I 0.22 (0.23) fm" 3 , for RHIC (b] 

with the MRSH (MRSD-') parton distributions. Computing the net baryon-to-entropy ratio 
by using s g = 3.6n g for a thermal boson gas gives initially, at r = 0.1 fm/c: (B — B)/ S g ~ 
2 • 10 -4 for LHC, and ~ 2 • 10~ 3 for RHIC. We conclude that at the future colliders we are still 
relatively far away from the extreme conditions of the Early Universe, where the inverse of 
the specific entropy is ~ 10 -9 [p7|] . For the LHC, assuming an instant thermalization of the 
gluon system at r = 0.1 fm/c, and an adiabatic evolution thereafter, the final entropy can 
be approximated by the initial entropy of gluons [|| [TT| . The non-perturbative mechanism 
for particle production will not increase the entropy much but does increase the net baryon 
number. If the non-perturbative contribution to the net baryon number production is as- 
sumed to be of the same order of magnitude as in the current Pb+Pb collisions at SPS, the 
final baryon-to-entropy ratio for the LHC will be ~ 10~ 3 . For the RHIC nuclear collisions, 
thermalization is most likely not as instantaneous, but following nevertheless the same line 
of arguments, and taking into account that the non-perturbative component becomes im- 
portant also for entropy production, we estimate ~ 10~ 2 for the final net baryon-to-entropy 
ratio. 



3 Minijets in the BFKL-approach 



Minijet production I have considered above is based on collinear factorization, where the per- 
turbative partonic cross sections are factorized at a momentum scale Q ~ pt from the parton 
distributions with nonperturbative input. Next, I will discuss an additional mechanism for 
minijet and transverse energy production, where factorization is not used. 

The small- £ rise in the structure function F 2 (x, Q 2 ) observed at HERA Q 2 > 1.5 GeV 2 
|22|j can be explained by the leading log(Q 2 ) DGLAP-evolution [2B| and also by the leading 
log(Q 2 )log(l/x) evolution g|. Also 

a power-like behaviour, F 2 ~ x , expected in the 
leading log(l/a;) BFKL-approach JHJ, does not contradict the data. In the following, let 
us assume that the small- x increase is entirely due to the BFKL-physics. Then, with this 
assumption, we will study what is the maximum transverse energy deposit in the central 
rapidity unit due to the minijets emitted from a BFKL-ladder in the LHC nuclear collisions. 
At RHIC the BFKL-minijets are not expected to contribute in any significant manner be- 
cause the BFKL-enhancement takes place only at x^ 0.01. Therefore, this latter part of my 
which is based on Ref. 



talk 
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will be relevant only for the LHC. 
It is instructive to start from a case of fully inclusive minijets with two tagging jets 

The (summed) 



separated by a wide rapidity g studied by Mueller and Navelet |31 



subprocess is also shown Fig. la, where the incoming partons have momentum fractions x a 
and Xb, the tagging jets rapidites y a and y^ (y a ^> yb) and transverse momenta k a x and k;/r, 
respectively. Between the tagging jets there are n gluons emitted, labeled by l...n. Thus 
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each final state is described by a Feynman graph with 2 incoming and n + 2 outgoing on-shell 
gluons. The colour singlet hard BFKL-pomeron ladder arises when these Feynman graphs 
are squared and summed. In the kinematic region we will be interested in, the rapidities are 
strongly ordered, y a ^> y\ 3> ... 3> y n 3> yb, but the transverse momenta are not, k Ti ~ k Tj -. 
Then only the transverse degrees of freedom of the momenta of the virtual legs become 
important. The tagging jets of Fig. la have transverse momenta at a perturbative scale, so 
that one may use collinear factorization to write the cross section down as: 



d 2 k aT d 2 ]s. bT dy a dy b ^ J H[ 4tt (2 

where only gluons are considered. The strong rapidity ordering simplifies the momentum 
fractions to x a ~ ^f eVa an d Xb ~ ^=e~ yb , and the parton densities factor out of the sum. 

For the process gg — > gg the leading contribution in the large s/t limit comes from 
the i-channel amplitude. In a physical gauge, this amplitude is also gauge invariant up to 
the subleading terms. The matrix element M a \_ nb consists then of the following building 
blocks: In the leading s/t approximation, in a physical gauge and with the strong rapidity 
ordering, each gluon can be regarded as emitted from an effective non-local Lipatov-vertex, 
where bremsstrahlung from initial and final legs and emission from the exchanged gluon 
are summed. These are described by the black blobs in Fig 1. Also the propagators are 
effective ones since they are exponentiated (Reggeized) after computing the leading virtual 
corrections to the t-channel gluon exchange. The effective propagators are drawn by thicker 
(vertical) lines in Fig. 1. Original references, detailed discussion and derivation of these 
concepts can be found in the useful lecture notes by Del Duca E2I . 



Next, we square each matrix element M aL .. n b, and due to the strong ordering in rapidities, 
colour singlet ladders with n + 2 rungs are formed. The colour algebra can be performed 
by summing (averaging) over the final (initial) state colours, and the polarization sums can 
be done. With help of e.g. Laplace-transformation (in y a — y^, the rapidity integrals can 
be disentangled. Finally, by summing over n, one obtains an iterative integral equation, 



the inhomogeneous BFKL-equation |30| £]Tj (see also |j32|), which describes an addition of 
one rung into the colour- singlet hard pomeron ladder. The BFKL-ladder is denoted by 
/(q T , k T , y a — y b ) in Fig la. The cross section (0) then becomes: 

da , 2\ / 2s 4A r c 2 «s 1 ,1 

-w, — -w, — , , = x a g{x a , fi ) x b g{x b , \i ) — — - 2/(q-r, k T , y a - Vb) tt, 
d 2 k aT d 2 k bT dy a dy b N 2 - 1 k A aT 



where qx = — k a x and kx = k^x- If all the virtual corrections and the real emissions are 
neglected, the ladder reduces into 2/(qx, kx, y) — > 5^(qx — kx), and the Born limit for the 
two jets separated by a large rapidity interval is recovered [ |3^| . 

Let us then study the case with tagging jets further by fixing one step of the ladder, as 
shown in Fig. lb. It is straightforward to sum the graphs with gluon emissions between the 
tagging jet a and the fixed minijet c, and, between the minijet c and the tagging jet b. This 
creates a ladder on each side of the fixed rung. Especially, we learn that a generic factor 
a s N c /(Tr 2 k^ T ), which includes the phase-space factor and contraction of the two Lipatov 
vertices associated with the step c, arises from fixing the the rung c. The cross section 
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,<Ut»J/« - Vc) 



,k T ,2/ c - ja) 



(a) (b) 

Figure 1: (a) The BFKL ladder in fully inclusive minijet production between the two tagging jets a and 
b ||, [33j. (b) Fixing one step of the ladder (c) creates a ladder on each side of the minijet c Q. 



becomes 



da 



d 2 k aT d 2 k bT d 2 k cT dy a dy b dy c 



x a g(x a , /j 2 ) x b g(x b , /1 2 



AN 2 a 2 a s N c 1 
N 2 -l~ 2 



-2- f c/ 2 qi T rf 2 q 2T 



d^(k cT - qi T + q 2T ) To n • ( 9 ) 



h 2 



Our goal is to study the leading BFKL minijet production mechanism which is ~ a s . 
As illustrated in Fig. 2, we therefore relax the requirement of having tagging jets. Then 
coupling of the pomeron ladder to the hadron becomes essentially non-perturbative and a 
form-factor, or, rather, a parton distribution, will be needed. Also, now that we do not 
require any tagging jets, we have to give up collinear factorization. We do not have any 
perturbative Born limit to compare with, either. Therefore, the best we can do is to adopt 



the procedure for deep inelastic scattering (DIS) in [35], where an addition of each rung into 
the pomeron ladder between the two hadrons or nuclei is expected to be described by the 
homogeneous BFKL equation for the unintegrated gluon density f(x, q\), 



— x- 



df(x, q\) a s N c 



dx 



71 



-Qt 



dq 



QlT 



f(x,q 2 T ) - f(x,q%) f(x,q% 



Qit\ 



+ 



It + 4 €t 



(10) 



Normalization for this scale-invariant equation is given by the gluon distributions 

"Q 2 dq\ 



xg(x, 



'ID 



determined from experimental input p5 |. 

By using the knowledge of the factor arising from fixing one rung of the BFKL-pomeron 
ladder, the inclusive minijet cross section from the BFKL-ladder can now be written down 
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Figure 2: Minijet production without the tagging jets requires an introduction of unintegrated parton 
distributions f(x,q^) Q- 



as 



dai ct 
d 2 p T dy 



K N ^—^ \ f d 2 q 1T d 2 q 2 T <5 (2) (Pt - qiT + Q2T\ 

TV 2 £4 J 



f(xi,qf) f(x 2 ,q%) 



Q 2 n 



<?2T 



(12) 



where p^ and y are the transverse momentum and the rapidity (in the hadron CMS) of the 
minijet. From momentum conservation and multi-Regge kinematics the momentum fractions 
become x X ( 2 ) ~ ^= e ±y . Due to the fact that in this case we do not have an "external" hard 
probe like the virtual photon with an associated quark box as in DIS, nor an on-shell Born 
cross section to relax into, we cannot determine the overall dimensionless normalization 
constant Kn exactly. However, we are able to fix the slope of the minijet pp-distribution, 
which will be sufficient for estimating the upper limit of transverse energy production from 
the BFKL-ladder. 

The minijet cross section of Eq. ([12]) is shown in Fig. 3 [14| . In the BFKL-computation we 



have used the unintegrated gluon densities compatible with the small-x rise in the set MRS- 
D-' |Jl7 |. For comparison, the more traditional (CFLTLO) minijet cross sections, discussed 
in the first half of the talk, are also shown with the MRSD-' parton distributions. The 
NLO jet analysis [|19|, [18| indicates that LO+NLO calculation with collinear factorization 
reproduces the measured jet cross sections well. Therefore, at 5 GeV, the BFKL-minijet 
contribution should be less than the collinearly factorized. We can thus argue that Kn< 1. 
The transverse energy production due to the minijets from the BFKL-ladders at \y\ < 0.5 



in AA collisions can now be estimated |T3 from 



dE* FKL 



dy 



T AA {h) 



,„BFKL 

Po 



dpTPT 



da iet 



dp T dy 



(13) 



The coherence of the BFKL ladder is broken when we fix a rung, and the cross section diverges 
at pt — > 0. A cut-off is, unfortunately, needed also in the BFKL case. The saturation of 
the CFLTLO-minijet cross section in the LHC Pb-Pb collisions (as considered in the first 
half of the talk), implies that the BFKL-cross section should not grow much larger than 
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BFKL: a,b; if iV =l, CFLTLO: 1,2,3 




1 2 3 5 10 

PT [OeV] 

Figure 3: The minijet cross sections at y = as functions of transverse momentum at \fs = 5.5 TeV 
from jl4j. (panel a) and y/s — 900 GeV (panel b) |Q. The curves 1, 2 and 3 are predicted in the approach of 
collinear factorization, leading twist and lowest order (CFLTLO) pQCD, and MRSD-' parton distributions. 
The curves a and b are the BFKL-minijets from Eq. ( |l2|) with and without running coupling. 



the curve 2 in Fig. |3|. Therefore, we do not trust the BFKL-computation with ~ 1 
below pt < Po FKL ~ 1 GeV. With these values, we find for central Pb-Pb collisions at the 
LHC, E% FKL = 3060 GeV with fixed a s = 0.2, and, 4940 GeV with (ad hoc) running a s (p T ). 
Comparing these numbers to the value 15330 GeV in Table lb for gluons, we see that the 
BFKL-contribution is at most a few 10% correction to the results in |4], |TTJ. On the other 
hand, one should perhaps compare the results at the same level of approximations, (only 
LO gluons, a s fixed) i.e. curves 3 and a in Fig. |^a. Then the two contributions become 
of similar magnitude. In this case, however, the po in the CFLTLO-computation should be 
lower than 2 GeV, and the BFKL contribution would again be the smaller one. 

We worked under the assumption that the BFKL-evolution is responsible for all the small- 
x rise at HERA, i.e. we studied the maximum contribution from the kinematical region 
relevant for the hard BFKL-pomeron. Since the HERA results can be explained by the 
leading \og(Q 2 ) and/or the leading log(Q 2 ) log(l/x) approximations, the leading log(l/x)- 
contribution is obviously not the dominant mechanism at the present values of x. Thus, 
my conclusion is that the BFKL-minijets certainly bridge the way towards softer physics at 
Pt < Po ~ 2 GeV, but the initial conditions relevant for the early QGP-formation in the LHC 
nuclear collisions are dominantly given by the minijets computed in collinear factorization. 
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